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Using numerical simulations, we investigate vortex configurations and pinning in superconductors 
with honeycomb and kagome pinning arrays. We find that a variety of novel vortex crystal states can 
be stabilized at integer and fractional matching field densities. The honeycomb and kagome pinning 
arrays produce considerably more pronounced commensuration peaks in the critical depinning force 
than triangular pinning arrays, and also cause additional peaks at noninteger matching fields where 
a portion of the vortices are located in the large interstitial regions of the pinning lattices. For the 
honeycomb pinning array, we find matching effects of equal strength at most fillings BjB^, = n/2 
for n > 2, where n is an integer, in agreement with recent experiments. For kagome pinning arrays, 
pronounced matching effects generally occur at B/B^, — n/3 for n > 3, while for triangular pinning 
arrays pronounced matching effects are observed only at integer fillings B /B^ — n. At the noninteger 
matching field peaks in the honeycomb and kagome pinning arrays, the interstitial vortices are 
arranged in dimer, trimer, and higher order n-mer states that have an overall orientational order. 
We call these n-mer states "vortex molecular crystals" and "vortex plastic crystals" since they 
are similar to the states recently observed in colloidal molecular crystal systems. We argue that 
the vortex molecular crystals have properties in common with certain spin systems such as Ising 
and n-state Potts models. We show that kagome and honeycomb pinning arrays can be useful for 
increasing the critical current above that of purely triangular pinning arrays. 

PACS numbers: 74.25.Qt 




I. INTRODUCTION 



There have been extensive studies on the static and 
dynamical properties of vortices in superconductors with 
periodic arrays of artificial pinning sites. These works 
focused on simple two-dimensional periodic pinning ar- 
rays such as square, triangular, and rectangular lat- 
tices, where the pinning sites consist of holes [J, H, H, 
blind holes [HI EH, or magnetic dots 
151 lip] . [T7L [l8j . For square and triangular pin- 
ning arrays, pronounced commensurability effects such 
as peaks or anomalies in the critical current appear at 
the magnetic field B — B^ where the number of vortices 
equals the number of pinning sites, as well as at higher 
fields B = nBtf,, where n is an integer. At these matching 
fields the vortex lattice forms ordered crystalline struc- 
tures of a type determined by the number of vortices 
that are captured at individual pinning sites. If more 
than one vortex can occupy each pinning site in the form 
of a multi-quanta vortex, then at each matching field the 
overall vortex lattice has the same symmetry as the pin- 
ning lattice but is composed of n-quanta vortices. If only 
a single vortex can occupy each pinning site, ordered vor- 
tex crystals still form at the matching fields and for n > 1 
some of the vortices are located in the interstitial regions 
between the pinning sites. Imaging experiments Q and 
simulations [T] for systems where at most one vortex can 
occupy each pin have shown that numerous kinds of inter- 
stitial vortex lattice structures can be stabilized, some of 
which have different symmetries than the pinning array. 
Interstitial vortex crystals also form above the pinning 



saturation field in samples with pins that can be occu- 
pied by multi-quanta vortices. When each pin has cap- 
tured as many vortices as possible, additional vortices sit 
in the interstitial regions, and the resulting vortex lattice 
structure is a composite of interstitial singly-quantized 
and pinned multiply-quantized vortices 0, M, Es . [20| . In 
the case where the pinning sites are blind holes, it is pos- 
sible for multiple vortices to occupy a single pinning site 
without merging into a multi-quanta vortex. Instead, the 
vortices retain their individual identities and form dimer, 
trimer, or ring- type states inside the pin [ToL [2ll| . 

In addition to the matching effects that appear at inte- 
ger fields, commensuration effects can also occur at non- 
matching fields or fractional fields I, US Eft Th ese 
noninteger matching effects are generally weaker than 
those observed at integer matching and are most promi- 
nent for fields B < B$. If multiple vortex quantization 
occurs above the first matching field, the sub-matching 
sequence that appears between B = and B — B^ is 
repeated between every integer matching field until the 
pinning sites are saturated. In contrast, if only one vor- 
tex is captured per pinning site the fractional matching 
effects above the first matching field are significantly re- 
duced or missing [8[ . Commensuration effects for vortices 
interacting with a periodic substrate have also recently 
been demonstrated for vortices in Bose-Einstein conden- 
sates where the pinning sites are created with an optical 
array HH^UgJ. 

The physics of vortices in periodic pinning arrays is 
similar to that of repulsively interacting colloids in trian- 
gular or square periodic trap arrays (2f| [26 . [27l . [2H 
and charged spheres on periodic substrates 3d||. In both 
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FIG. 1: The pinning site locations (circles) for periodic ar- 
rangements of sites forming (a) triangular, (b) kagome, and 
(c) honeycomb arrays. Here — O.470o/A 2 . 



these cases it is possible to have localized traps which 
capture only a single colloid or sphere while the remain- 
ing particles sit in the interstitial regions [29L l3fjl ] , similar 
to the situation for vortex pinning arrays. The interstitial 
particles are more mobile than the pinned particles, and 
the particle trajectories resemble those seen in computer 
simulations of vortices in similar geometries 
the colloid system it is also possible for multiple colloids 
to be captured by a single trap (2f| [26|, [U HI] where they 
form dimer, trimer, or higher order n-mer states. The ri- 
mers may be orientationally ordered if the interaction 
between colloids in neighboring traps is strong enough, 
and the system can show multistage melting transitions 
in which the initial stage of melting occurs when the ori- 
entational ordering of the n-mers is lost, followed by a 
transition in which the n-mers break apart. The ordered 
n-mer states can be mapped onto various spin systems 
such as Ising and Potts models [13, • Orientationally 
ordered dimer and trimer states have also been proposed 
to occur for vortices in superconductors with blind hole 
arrays where the vortices retain single quantization in 
an individual hole [8(. The orientational ordering of the 
n-mer states in the colloidal and vortex systems arises 
due to quadrupole or higher order pole moment inter- 
actions between the n-mers in adjacent traps [27]. The 
anisotropic nature of the quadrupole and higher order 
pole moments creates preferred directions for the align- 
ment of the n-mers that minimize the pole moment en- 
ergy. 

A limited number of studies have treated periodic sub- 
strates other than square and triangular pinning arrays, 
such as honeycomb or kagome pinning arrays. These are 
simply triangular lattices with a fraction of the pinning 
sites removed. Figure [IJa) shows a triangular pinning 
lattice. Removing every other pinning site from every 
other row of the triangular lattice, which eliminates 1/4 
of the pinning sites, produces the kagome array illus- 
trated in Fig. HUb) . A honeycomb array in which 1 /3 of 
the pinning sites in the triangular array have been elim- 
inated is shown in Fig. [He). Here every third pinning 
site is removed from every row of the triangular array. 
An experimental study of pinning phenomena in kagome 
arrays using superconductors with magnetic dot arrays 
produced evidence for pronounced commensurability ef- 



fects at noninteger matching fields [l3(. These results 
are, however, difficult to interpret since magnetic dot ar- 
rays can induce the formation of vortices or antivortices 
in addition to the vortices created by the external field 
(l7| . Numerical work on kagome pinning arrays [3l| only 
treated thermal melting of vortices at a single field of 
B/B^ = 2.0 where a two step melting transition was 
shown to occur. Here, the multiple interstitial vortices 
located at the larger interstitial sites where a pin has been 
removed undergo local melting at a temperature below 
that at which the entire vortex lattice melts. In Ref. [3l| . 
neither long range vortex configurations nor critical cur- 
rents were analyzed, so the nature of commensurability 
effects for varied vortex densities in kagome pinning ar- 
rays is not known. 

Only a single experimental study (33j has been per- 
formed on honeycomb pinning arrays to our knowledge. 
In Ref. [331 ]. several unusual features were observed, in- 
cluding pronounced matching effects of equal magnitude 
at magnetic fields corresponding to both integer and half 
integer multiples of the matching field for vortex densities 
up to to the fifth matching field B j = 5. This result is 
in contrast to the response of square or triangular pinning 
arrays, where commensuration effects are much weaker 
at nonmatching fields than at integer matching fields. In 
Ref. HU , it was also observed that for fields greater than 
the second matching field, B/B^ > 2, the commensura- 
tion effects at half-integer matching fields become more 
prominent than those at integer matching fields. This 
effect can be understood by considering that since the 
honeycomb pinning array is simply a triangular pinning 
array that has been diluted by 1/3, the field B/B^ = 1.5 
for the honeycomb array would correspond to an integer 
matching field B/B^ = 1 in a triangular array of the 
same density. As a result, the overall vortex lattice is 
triangular at B / B^ = 1.5 in a honeycomb pinning array. 
Similarly, a field of BjB^ = 1/2 in the honeycomb pin- 
ning array would correspond to a field of B / B^ = 1/3 in 
the equivalent triangular pinning array, which is known 
to produce a peak in the critical current @. These re- 
sults suggest that honeycomb pinning arrays may allow 
for a variety of new vortex structures to be stabilized at 
noninteger fillings. 

In this work we present the first extensive study of 
vortex pinning and dynamics in honeycomb and kagome 
pinning arrays using numerical simulations. Section II 
contains a description of the simulation method. We 
show the vortex configurations and ordering for honey- 
comb pinning arrays in Section III and illustrate the for- 
mation of vortex molecular crystals, which are named in 
analogy with molecular crystals. Section IV gives the cor- 
responding description for vortices on a kagome pinning 
lattice. The melting of these vortex configurations and 
the creation of vortex plastic crystal states for both types 
of pinning lattices is studied in Section V. We construct 
phase diagrams for the vortex molecular crystals at the 
dimer and trimer fillings of the honeycomb pinning lat- 
tice as a function of temperature and pinning strength in 
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Section VI. The effect of the strength of the pinning sites 
is explored in further detail in Section VII. The paper 
closes with a discussion in Section VIII and a conclusion 
in Section IX. 



II. SIMULATION 

We perform two-dimensional simulations of supercon- 
ducting vortices in honeycomb and kagome pinning ar- 
rays using a computational procedure similar to that pre- 
viously employed for studies of vortices in square and tri- 
angular pinning arrays @, H, HH • The system of size 
L x x L y has periodic boundary conditions in the x and 
y directions and contains N p pinning sites and N v vor- 
tices. The dynamics of a vortex i located at position R; 
is determined by the following overdamped equation of 
motion: 



n- 



~dt~ 



(1) 



Here F™ is the repulsive vortex-vortex interaction force, 
is the force from the pinning sites, Ff is the force 
from an external drive, and Ff is the random force 
from thermal fluctuations. The damping constant r/ = 
cf>od/2ir£ 2 pm where d is the thickness of the supercon- 
ducting sample, <fio = h/2e is the flux quantum, £ is the 
superconducting coherence length, and p^ is the normal 
state resistivity of the material [34j . We measure length 
in units of the London penetration depth A and for most 
of the results presented here the system size is 24A x 24A. 
The explicit vortex-vortex interaction force is 



R 



> J ! 



(2) 



where K\ is the modified Bessel function, fo = 
<j> 2 /(2np \ 3 ), Rij = |Ri-Rj|, and fty = iR, R, i //,.,. 
A short range cutoff of 0.1A is applied to the vortex- 
vortex interaction force to avoid divergence; however, at 
the densities considered here, vortices do not approach 
each other this closely. The interaction is also cut off 
beyond 6A for computational efficiency since the vortex- 
vortex forces beyond this distance are negligible Q • Time 
is measured in units of r = T]/ fo- As an example of phys- 
ical units, for a NbSe2 crystal of thickness d = 0.1 mm 
with rj = 2.36 x lCT 11 Ns/m, f = 6.78 x 1CT 5 N/m and 
t = 0.35 ps. 

The pinning sites are modeled as N p attractive 
parabolic traps with radius R p = 0.3A and strength f p , 
so that 



R n — R 



(p)' 

ik 



k=l 



Ru 



R 



(v) 



(3) 



Here R^ is the location of pinning site k, R^ 



|Ri 



(Ri - B^ } )/r\1\ and 6 is the Heaviside 



step function. The pinning radius is chosen to be small 
enough that only one vortex can be captured per pinning 
site. The pinning sites are arranged in a triangular lattice 
with lattice constant ao, as illustrated in Fig. [3a), and 
then 1/4 or 1/3 of the pinning sites are removed to create 
a kagome or honeycomb array, as in Fig. [T]Jb,c). The 
pinning density is n p — N p / (L x L y ). The matching fields 
for the honeycomb and kagome arrays are B^f = n p and 
B^ = n p , respectively, while the matching field of the 
equivalent triangular array is -B^, such that B? /B^ = 
2/3 and B* /B = 3/4. 

The driving force is assumed to arise from the appli- 
cation of an external current which induces a Lorentz 
force on the vortices that is perpendicular to the current. 
All vortices experience an equal driving force F d = F d x 
in the x direction, corresponding to the horizontal axis 
of Fig. [TJ The thermal force Ff is modeled as ran- 
dom Langevin kicks with the properties (F^"} = and 
(Ff (t)Fj(t'))=2 V kBT5(t-t>)5 ij . 

The initial vortex configurations are obtained by sim- 
ulated annealing. Our procedure for this study was to 
start from an initial temperature of F T = 3 and decrease 
the temperature to F T — in increments of SF T = 0.002 
while spending 5000 simulation time steps at each in- 
crement, so that the total annealing time is 7.5 x 10 6 
simulation time steps. We note that an overly rapid an- 
nealing rate can cause the system to be trapped in a 
metastable state, which prevents the vortices from or- 
dering even at integer matching fields. To check our an- 
nealing rate, we tested slower rates and found that the 
resulting vortex configurations were unchanged. Once 
the vortex positions have been initialized and the tem- 
perature has been set to zero, we determine the velocity- 
force curve relations and the critical depinning force f c by 
slowly increasing the external drive F d . We measure the 
average vortex velocity response V x = N' 1 ^^ v, • x), 
where Vj = dUi/dt. The resulting velocity- force curve 
would correspond to a voltage-current curve in experi- 
ment. The depinning force is determined by applying a 
cut-off threshold of V x = 0.01 to the average velocity. We 
find that this cutoff is sufficiently low that the critical de- 
pinning force f c vs B curves are not strongly sensitive to 
the choice of cutoff. 



III. VORTEX PINNING AND ORDERING IN 
HONEYCOMB ARRAYS 

A. Commensurability Peaks at Integer and Half 
Integer Fillings 

We measure the critical depinning force f c for vor- 
tices on a honeycomb pinning array with f p — 0.5/o and 
B0 = 0.47</>o/A 2 as a function of vortex density. The 
results are plotted in Fig. ^a.) as f c /f P versus B/B? 
where B? = 0.313<^o/A 2 is the field at which the num- 
ber of vortices equals the number of pinning sites for the 
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FIG. 3: (Color online) f c /f P vs BjB$ for the honeycomb 
pinning array (black heavy line) and the triangular pinning 
array (red light line) in Fig. l(a,c) with f p = 0.5/o and B$, = 
0.47<^o/A 2 . 
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FIG. 2: (a) The depinning force f c /f p vs B/Bf for the 
honeycomb pinning array illustrated in Fig. die) with f p — 
0.5/o, B<j, = O.470 o /A 2 , and Bf = O.3130 o /A 2 . (b) A blow-up 
from (a) of f c /f p vs B/Bf in the region B/Bf > 2. 



honeycomb pinning array. We find peaks in f c at the 
integer matching fields B/B^ = 1, 2, 3, and 4. Addi- 
tionally, there are pronounced peaks in f c at the half- 
matching fields B/B% = 0.5, 1.5, 2.5, 3.5, and 4.5. Note 
that the f c peaks at the noninteger matching fields of 
B/B^ = 3.5 and 4.5 are significantly larger than the 
peaks at the integer matching fields of B/B^ — 3, 4 and 
5, as highlighted in Fig[2jb). There are sub-matching 
commensuration effects for 0.5 < B/B^ < 1.0 which are 
in general weaker than the sub-matching commensura- 
tion effects for B/B^ = and 1.0. A similar trend was 
observed in earlier works for square and triangular pin- 
ning arrays [H, Q. The numerical simulations are time 
consuming and permit us to perform only finite field in- 
crements, so that the weaker submatching and higher 
order fractional matching fields are difficult to observe 
numerically. In this paper we focus mainly on the pro- 
nounced matching fields. 

The behavior of the commensuration effects for the 
honeycomb pinning lattice in Fig. [2] is very similar to 
the experimental results observed for honeycomb arrays 
in Ref. [33j . In the experiments, strong commensura- 
tion effects appeared at B/B^ = 1/2, 1, 1.5, and 2, in 
agreement with our results. The experiments also show 



that the commensuration effects at B/B^ = 3 and 4 
are very weak or absent while those at B/B^ = 3.5 and 
4.5 are very pronounced. This is also in agreement with 
our result as seen in Fig. [2jb). One difference between 
the experimental results and our work is that we find a 
very pronounced peak in f c at at B/B^ = 2.5, while 
the experiments produced only a weak peak at the same 
field. In Section V we show that this may be due to 
the fact that the commensurate vortex configuration at 
B/Bf = 2.5 is unstable under thermal fluctuations. 



Figure [3] shows a comparison between the behavior of 
f c as a function of B for the honeycomb pinning array 
(heavy line) and a triangular pinning array with the same 
value of B^ = 0.47</>o/A 2 (light line). The depinning force 
f c for the triangular pinning array exhibits pronounced 
peaks at B /B^ = 1/3,2, and 3, and only very weak peaks 
at noninteger matching fields for B/B^ > 1. The peak 
at the second matching field B/B^, = 2 for the trian- 
gular pinning array is significantly smaller than the one 
at the third matching field B/B^ = 3. This effect has 
been observed in previous studies and results from the 
fact that the vortex lattice forms a honeycomb structure 
at the second matching field while at the third match- 
ing field the overall symmetry of the vortex lattice is 
triangular [7| . Figure [3] also indicates that the peak at 
BjB^ = 1.0 for the triangular array coincides with the 
peak at B/B^ = 1.5 for the honeycomb array, while a 
peak at B/B^ = 1/3 in the triangular array coincides 
with a peak at B /B? = 1/2 in the honeycomb array. 
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FIG. 4: Left panels: Vortex positions (black dots) and pin- 
ning site positions (open circles) obtained for the honeycomb 
array in Fig. Die) an d Fig. [2] Right panels: Vortex posi- 



tions only. (a,b) B/Bf 
B/B% = 1.5. 



1/2. (c,d) B/B" = 1. (e,f) 



B. Vortex Configurations for B < B^, 

In order to explain the origin of the pronounced com- 
mensurability peaks in f c for the honeycomb pinning ar- 
ray at B/B^ = n/2 and the fact that some peaks are 
more prominent than others, we analyze the vortex con- 
figurations at these fields. In the left panels of Fig. [4] we 
illustrate the vortex and pinning site locations for three 
matching fields in the honeycomb array from Fig. [21 while 
the right panels of Fig. 0] show only the vortex positions 
at the same fields in order to emphasize the vortex lattice 
structure. 

Figure Ufa, b) indicates that at B/B? = 1/2 the over- 
all vortex lattice is triangular. The vortex configuration 
at half filling for the honeycomb lattice is the same as 
that of vortices at unit filling in an undiluted triangular 
lattice with matching field £?? /2 that is rotated by 90° 



relative to the triangular lattice in Fig.QJa). The config- 
uration is also identical to the vortex arrangement found 
for 1/3 filling of a triangular lattice with matching field 
B<f, |8j. Since the vortex lattice structure is triangular, 
the vortex-vortex interactions cancel and the depinning 
force is determined only by the maximum force of the 
pinning sites such that f c / f p = 1.0. This is what we find 
at B/B^ = 1/2, as seen in Fig.fS] For fields just above or 
below half filling, the vortex lattice retains the same tri- 
angular ordering shown in Fig.|3]Ja,b) but contains weakly 
pinned vacancies or interstitials which reduce the value 
of f c . In the case of a triangular pinning array, there is 
no commensurate peak at B/B^ = 1/2 when all of the 
vortices are forced to occupy pinning sites since the sys- 
tem is geometrically frustrated, resulting in a strongly 
defected vortex lattice [8j. If f p is very weak an ordered 
vortex lattice can form at this field when the elastic forces 
of the vortex lattice overcome the pinning force, allowing 
half of the vortices to shift out of the pinning sites and 
creating a partially pinned or floating triangular vortex 
lattice [8|, l3fj |. In our simulations the pinning strength 
f p = 0.5/o is well above this limit so that floating vortex 
configurations do not occur; we return to this point in 
Sections VI and VII. We note that for very low applied 
fields B/B,/, <C 1, the depinning force is dominated by 
single vortex pinning and thus f c /fp ~ 1- 

The vortex configuration for the first matching field 
in the honeycomb array, B/B^ — 1, is illustrated in 
Fig. Etc,d). The vortex lattice has the same highly sym- 
metric structure as the pinning lattice and as a result the 
vortex-vortex interactions cancel, giving f c /f P — 1-0 as 
shown in Fig. [2] For fields above the first matching field, 
the additional N v — N p vortices are located in the inter- 
stitial regions and are pinned not by the pinning sites 
but by the interactions with the vortices trapped at the 
pinning sites. In general this interstitial pinning is weak, 
so fc/fp drops significantly for B/B^ > 1.0 as seen in 

Fig.m 

At B/B^ = 1.5, where a prominent peak in the de- 
pinning force appears in Fig. [21 the overall vortex lattice 
is triangular as shown in Fig. 0|e,f). In this case the 
interstitial sites that were produced when the triangular 
pinning array was diluted to form the honeycomb pinning 
array each capture one vortex. Figure E^b) shows that 
the prominent peak in f c /f P at B / B^ — 1 for the triangu- 
lar pinning array coincides with the peak at B/B^ = 1.5 
in the honeycomb pinning array. Although the symme- 
try of the vortex lattice is triangular in each case, at 
B/B^ = 1.5 for the honeycomb pinning array the depin- 
ning force is determined not by f p but by the caging force 
on the vortices in the interstitial regions. As a result, the 
depinning force for the honeycomb array at this field is 
lower than the depinning force for the triangular array 
at B/B^ = 1.0 filling, as seen in Fig. [3fb). For fields 
slightly below or above B/B^ = 1.5 in the honeycomb 
pinning array, interstitials or vacancies appear in the tri- 
angular vortex lattice at the locations of the missing pins 
and cause a reduction in f c / f p . 
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FIG. 5: Left panels: Vortex positions (black dots) and pin- 
ning site positions (open circles) in a portion of the sample 
for the honeycomb pinning array from Fig. [He) and Fig. [2] 
Right panels: Vortex positions only. (a,b) B/B^ = 2.0. (c,d) 
B/Bf = 2.5. (e,f) B/Bf = 3.0. 



C. Vortex Configurations for B > B$ 

As B increases above B/B^ =1.5, the additional vor- 
tices sit in the large interstitial regions at the center of 
each honeycomb plaquette. In general for 1.5 < B/B^ < 
2.0 we find that the plaquette centers capture at most 
two vortices rather than three. The two vortices cannot 
both sit at the center of the interstitial site, so instead 
they form an interstitial dimer state, as illustrated in 
Fig. [Ua,b) at B/B^ = 2.0. The interstitial dimers have 
an additional orientational ordering which is highlighted 
in Fig. OJa). The dimers can be described as having a 
director field which is oriented at a 30° angle with the 
x axis in Fig. EJa,b). Within the center interstitial re- 
gion of each honeycomb plaquette, the six surrounding 
pinned vortices create a sixfold modulated symmetric po- 



tential and each vortex forming the interstitial dimer sits 
in one of the minima of this potential. As a result, at 
zero temperature each dimer can be oriented in one of 
three degenerate directions. If we change the initial con- 
ditions during the simulated annealing process slightly, 
we find the same configuration shown in Fig. E)^a,b) with 
one-third probability. With equal probabilities we obtain 
configurations in which all the dimers are either oriented 
at 90° or —30° angles with the x axis. If the dimers in 
neighboring plaquettes did not interact with each other, 
we would expect to find a random distribution of dimer 
orientations among the three degenerate directions in a 
given configuration. The orientational ordering of the 
dimers in our system indicates that dimers in neighboring 
plaquettes do interact with each other, and that this in- 
teraction gives the dimers a ferromagnetic-like alignment. 
Unlike an Ising model in zero magnetic field which has 
two possible spin orientations, this system has three pos- 
sible orientations for the dimers and is thus more closely 
related to a three-state Potts model. The ordering of the 
dimers is very similar to the ordering found in the re- 
cently studied model of colloidal dimers on a triangular 
lattice, which has been shown to map to the three-state 
Potts model [H]. The theoretical work in Ref. [H] in- 
dicates that the dimers lose their orientational ordering 
as a function of temperature through either a contin- 
uum melting transition or a first order transition, de- 
pending on the system parameters. The colloidal dimers 
in neighboring plaquettes have been shown to interact 
through an effective quadrupole moment with an addi- 
tional screening term [27j | ; higher order n-mers were also 
considered which interact through higher order pole mo- 
ments. These types of interactions are anisotropic and 
thus the pole moment energy can be minimized when the 
n-mers form an orientationally ordered state. Neighbor- 
ing dimers may be oriented parallel or perpendicular to 
each other depending on the pinning geometry [25] . The 
colloids interact via a repulsive screened Yukawa poten- 
tial; since this is similar to the vortex interaction of Eq. 
2, the same type of multipole interactions between vor- 
tices in neighboring plaquettes should emerge as in the 
colloid system. 

For fields 2.0 < B/B^ < 2.5, the additional vortices 
again occupy the interstitial sites where they form trimer 
states. For this range of fields we do not find any inter- 
stitial sites that have captured four vortices. Figure [HJc) 
shows the vortex configuration at B/B^ = 2.5 where 
each pinning site captures one vortex and each center in- 
terstitial region captures three vortices. The trimers are 
equilateral triangles and each vortex in the trimer is lo- 
cated at one of the sixfold potential minima created by 
the six surrounding pinned vortices. In the same man- 
ner as the dimer states, the trimers are orientationally 
ordered and all align in one of the two possible degener- 
ate orientations, indicating that neighboring trimers have 
an interaction with ferromagnetic character. The overall 
vortex lattice structure at B/B^ = 2.5 is very intricate, 
as indicated in Fig.[5£d). It can be viewed as triangles of 
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vortices each surrounded by three pentagons, indicating 
that a large fraction of the system has fivefold ordering 
coexisting with true long range order. 

Figure EJe) illustrates the vortex configurations and 
pinning site locations at B/B^ = 3.0. At this filling a 
weak peak in f c /f P appears, as shown in Fig. [2] Here each 
pinning site captures one vortex and the center intersti- 
tial regions capture four vortices in a quadrimer state. 
Unlike the dimers and trimers, we find no orientational 
ordering of the quadrimers. The fourfold symmetry of the 
quadrimer state cannot match the sixfold modulation of 
the potential at the center of the honeycomb plaquette. 
If the four vortices in the quadrimer sit at minima of this 
potential, then at least two vortices must occupy adja- 
cent minima, which is energetically unfavorable due to 
the vortex-vortex interactions. Instead, the quadrimers 
form distorted square structures that are not commen- 
surate with the underlying sixfold modulated potential. 
Since the overall orientational order of the quadrimers is 
absent at the B/B^ = 3.0 filling, a strong commensu- 
rate peak in f c / f p does not occur. A weak peak does 
appear at this filling, as seen in Fig. ^ due to the fact 
that each interstitial region captures exactly four vortices 
at B/B^ = 3.0. Just above or below this field, vacancies 
or interstitials in the form of threefold or fivefold occu- 
pied interstitial sites occur which are less strongly pinned 
and cause a reduction in f c . 

In Fig. [5Ja,b) we show the vortex configurations at 
B/B^ = 3.5 where a strong peak in f c /f P appears in 
Fig. [DJb). At this field, the system deviates from the 
pattern which we observed at lower fields for n = 1 to 
4 of forming symmetrical n-mers in the interstitial sites 
at the center of the honeycomb plaquettes. Instead of 
forming a pentamer state at B/B^ = 3.5, the intersti- 
tial vortices are arranged with four vortices captured as 
a rectangular quadrimer in the center of each interstitial 
plaquette and a fifth vortex between two adjacent pinned 
vortices. The vortex lattice has long range order with all 
the center quadrimers aligned in the same direction, un- 
like the B/Bg = 3.0 filling of Fig.Efof). The ordering is 
possible because the fifth interstitial vortex which sits at 
the boundary of the honeycomb plaquette breaks the six- 
fold symmetry of the potential minima inside the center 
of the plaquette and replaces it with a twofold symme- 
try which can be matched by the remaining quadrimer 
of interstitial vortices. 

At B/B^ = 4.0, a small commensuration peak ap- 
pears in Fig. 2(b). The corresponding vortex configu- 
ration is illustrated in Fig. [Hc,d). The vortex lattice 
has long range order. Each interstitial site captures six 
vortices in the form of an inverted triangular structure 
with three interstitial vortices at the top, two in the 
middle and one at the bottom. The triangles have a 
twofold degenerate ordering, and we have also observed 
the other possible ordering in which each triangle points 
upward instead of downward along the y direction. Al- 
though the vortex lattice structure has long range order 
at B/Bf = 4.0, there is not a large increase in f c /f P at 
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FIG. 6: Left panels: Vortex positions (black dots) and pin- 
ning site positions (open circles) in a portion of the sample for 
the honeycomb array from Fig. U^c) and Fig. [2] Right panels: 
Vortex positions only. (a,b) B/Bf = 3.5. (c,d) B/B" = 4.0. 
(e,f) B/B$ = 4.5. 



this field. This may be due to the fact that the trian- 
gular configuration of interstitial vortices is unstable to 
fluctuation effects, as we will describe in Section V. 

In Fig. EIe,f) we illustrate the vortex configuration at 
B/B^ = 4.5 where a pronounced peak in the depinning 
force occurs in Fig. [2h). In this case the seven inter- 
stitial vortices captured in each interstitial site are ar- 
ranged with one vortex in the center surrounded by six 
interstitial vortices sitting in the sixfold potential min- 
ima created by the six neighboring pinned vortices. The 
overall vortex lattice has a triangular ordering as seen 
in Fig. [6{f), and the vortex configuration is the same as 
that which would result at B/B^ = 3.0 for a triangular 
pinning array that is rotated by 90° relative to the trian- 
gular lattice in Fig. [TJa) , where the center interstitial vor- 
tex would be located in a pinning site. For B/Bf = 5.0 



(not shown) we find a state without orientational order 
and there is no particular peak in f c /f p at this field. We 
note that for very high applied fields B/Bf 3> 1, the 
dcpinning response is dominated by shearing motion of 
the interstitial vortices. 

As we have demonstrated, the honeycomb array can 
stabilize various types of interstitial vortex n-mer states 
at the center of each honeycomb plaquctte. These n-mers 
have a tendency to align in the same direction, indicat- 
ing effective ordering of a ferromagnetic nature which is 
similar to the ordering observed for colloidal n-mer states 
on periodic substrates. At fields B/Bf = n/2 where an 
ordered vortex crystal forms we observe peaks in the de- 
pinning force. As we describe in further detail in Sec. V, 
we find that at finite and increasing temperature, a melt- 
ing transition can occur in which the n-mers lose their 
orientational ordering but remain trapped in the center 
interstitial regions of each honeycomb plaquette, similar 
to the vortex dimer states in kagome arrays [3l[ and the 
colloidal n-mer molten states [2j| [2(| [28| . This suggests 
that many of the spin model analogies developed to de- 
scribe orientational ordering of colloidal n-mer states can 
be applied to the vortex n-mers as well. We term the ori- 
cntationally ordered vortex n-mer states "vortex molec- 
ular crystals" in analogy with molecular crystals, which 
have translational order along with an additional align- 
ment of the molecules. At higher temperatures where 
the molecules lose their orientational ordering but re- 
main in translationally ordered lattice positions, the sys- 
tem is referred to as a plastic crystal. Thus, by analogy, 
states such as that in Fig. [5^e) would be a vortex plas- 
tic crystal. The high temperature states (described in 
Sec. V) in which the orientational ordering has melted 
at B/Bf = 2, 2.5, and 3.5 would also be plastic vortex 
crystals. 



IV. VORTEX STATES AND 
COMMENSURABILITY IN KAGOME ARRAYS 



A. Commensurability Peaks at Integer and 
B/Bf = n/3 Fillings 

In Fig. [TJa) we plot f c / f P obtained for vortices in- 
teracting with the kagome pinning array illustrated in 
Fig. Htb) with f p = 0.5/o, B^ = O.470o/A 2 , and Bf = 
O.35250 o /A 2 , where Bf is the field at which the num- 
ber of vortices equals the number of pinning sites in the 
kagome array. Unlike the honeycomb array, there is no 
strong peak in f c /f P at B/Bf — 1/2 or at any of the 
submatching fields B/Bf < 1 for the kagome array, but 
clear commensurability peaks occur at the integer match- 
ing fields B/Bf — 1, 2, 3, and 4. In addition to the inte- 
ger peaks there are a series of peaks at Bj Bf — n/3 for 
n > 3. The strongest of these peaks fall at B/Bf =4/3, 
5/3, 7/3, and 10/3, as shown in Fig. Efb). The exis- 
tence of clear fractional matching effects at fields above 
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FIG. 7: (a) The depinning force f c /f p vs B/Bf for the 
kagome pinning array illustrated in Fig.rjjb) with f p = 0.5/o, 
B<j> = 0.47<?!>o/A 2 , and Bf = 0.3525<^o/A 2 . (b) A blowup from 
(a) of f c /f p vs B/Bf in the region B/Bf > 2.0. 

the first matching field is similar to the behavior seen 
in the honeycomb pinning array (Fig. [2]) and distinct 
from the behavior of a triangular pinning array (Fig. [3]), 
where no strong fractional matching effects appear above 
B I Bfj, = 1. This indicates that the pinning behavior of 
the kagome and honeycomb pinning arrays is very simi- 
lar and suggests that similar types of orientationally or- 
dered vortex molecular crystal states are occurring in the 
kagome array as in the honeycomb array. 

B. Vortex Configurations for B < B^, 

In Fig. [8{a,b) we show the vortex configuration for the 
kagome pinning array from Fig. [7]at B/Bf = 1/3. Here 
the vortex lattice does not order and in general we do not 
observe any particularly ordered vortex lattices for the 
sub-matching fields B/Bf < 1 for the kagome pinning 
arrays, in contrast to the triangular and honeycomb pin- 
ning arrays. In the honeycomb pinning array, Fig. [2]Ja) 
showed that there is a peak in f c /f P at B/Bf = 1/3 
followed by a drop in f c /f P - For the kagome pinning, 
Fig. [7] indicates that although there is no peak in f c /f p 
at B/Bf = 1/3, there is still a drop in f c /f p at this field. 
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FIG. 8: Left panels: Vortex positions (black dots) and pin- 
ning site locations (open circles) obtained for the kagome 
array in Fig. [TJb) and Fig. [7] Right panels: Vortex posi- 
tions only. (a,b) B/Bf = 1/3. (c,d) B/Bf = 1.0. (e,f) 
B/Bf = 4/3. 



In each case, the drop in f c /fp occurs due to a change in 
the nature of the vortex- vortex interactions. For fields at 
and below 1/3 filling, the spacing between adjacent vor- 
tices is at least 2ao and the vortex-vortex interactions are 
minimal. For fields above 1/3 filling, in order for all of 
the vortices to occupy pinning sites, some of the vortices 
must be located at a distance of a single lattice constant 
clq from another vortex, while the spacing between other 
nearest-neighboring vortices remains the larger distance 
2aQ. A vortex which has some nearest neighbors at a 
distance ao and other nearest neighbors at a distance 
2ao experiences an asymmetric vortex- vortex interaction 
force. This asymmetry causes the vortex to depin at a 
significantly lower driving force, producing the drop in 
fc/fp above 1/3 filling. 

At B/Bf = 1.0 where each pinning site captures one 
vortex, the net vortex symmetry is that of a kagome lat- 
tice, as indicated in Fig. |5fc,d). At B/B¥ — 4/3, shown 



in Fig. EKe,f), each center interstitial site of the kagome 
plaquettes captures one vortex so that the overall vor- 
tex configuration is triangular. This vortex configuration 
is the same as that for the triangular pinning lattice in 
Fig. QJa) at B/B^ = 1, where the vortices at intersti- 
tial locations in the kagome lattice would sit in pinning 
sites in the triangular pinning lattice. The B/Bf = 4/3 
filling also corresponds to the B/Bf = 1.5 filling of the 
honeycomb array where each interstitial site captures one 
vortex, and to the "first" matching field of Ref. [3l| . For 
the kagome lattice at Bj Bf <C 1 , single vortex depinning 
dominates the response of the system and f c / f p « 1. 



C. Vortex Configurations for B > B^ 

Figure OJa,b) illustrates the vortex configuration at 
B/Bf — 5/3 in the kagome pinning array. At this filling, 
each interstitial site at the center of the kagome plaque- 
ttes captures two vortices which form a dimerized state 
as highlighted in Fig. [Ha). Unlike the dimerized state 
in the honeycomb pinning array where the dimers are all 
aligned in the same direction, in the kagome array the 
dimers are tilted in opposite directions from one row to 
the next. This type of dimer structure is referred to as 
a herringbone state and has been observed for colloidal 
dimers on triangular arrays [25|, [28| as well as for the de- 
position of molecular dimers on triangular substrates [HJ 
and in three-state Potts models [36| . 

At B/Bf = 2.0 each central interstitial site captures 
three vortices which form a trimer state, as shown in 
Fig. [9jc,d). The trimer orientation is twofold degenerate 
as in the honeycomb array, and the trimers are oriented 
pointing either up or down. Fig. [5]Jc) shows that the up- 
down ordering is neither random nor uniform but occurs 
in domains reminiscent of an Ising model in zero field 
where domain walls have formed. If the domain wall 
energy is low, the relaxation time to reach a state where 
one of the phases dominates can be exceedingly long, and 
here we find that it is beyond our computational time 
scale. In spite of the lack of true long range order, a 
small peak still appears in f c / f p at this filling, as seen in 

Fig. raw. 

In Fig. M,e,i) at B/Bf = 7/3, the large interstitial 
sites in the kagome plaquette centers again capture three 
vortices and the additional fourth interstitial vortex is lo- 
cated in the small interstitial space between three pinning 
sites. Only half of the small interstitial regions capture a 
vortex. The trimers in the large interstitial sites are all 
aligned in the same direction, in contrast to the mixture 
of alignments found at B/Bf = 2.0 in Fig. E£c,d). The 
presence of the fourth interstitial vortex in the small in- 
terstitial sites breaks the sixfold symmetry of the large 
interstitial region and produces only a single low-energy 
alignment direction for the interstitial trimer. A grain 
boundary would require a shift in the position of the 
vortices in the small interstitial sites as well as a change 
in the orientation of the trimers. Such a grain bound- 
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FIG. 9: Left panels: Vortex positions (black dots) and pin- 
ning site locations (open circles) in a portion of the sample 
for the kagome array from Fig. [7J Right panels: Vortex po- 
sitions only. (a,b) B/Bf 
B/Bf = 7/3. 



FIG. 10: Left panels: Vortex positions (black dots) and pin- 
ning site locations (open circles) in a portion of the sample 
for the kagome array from Fig. [7] Right panels: Vortex po- 
5/3. (c,d) B/Bf = 2.0. (e,f) sit ions only. (a,b) B/Bf = 8/3. (c,d) B/B$ = 3. (e,f) 

B/Bf = 10/3. 



ary has a high energy and is not stable in our simulation. 
Since there are no grain boundaries present, this state has 
long range orientational order, and as a result the peak 
in fc/fp at B/B^ = 7/3 shown in Fig. [7] is higher than 
the peaks at B/B^ = 2.0 and B/B£ = 8/3 where the 
vortices do not form completely oricntationally ordered 
states. 

A related state appears at B/B^ = 8/3, shown in 
Fig. fTOTa.b). where the large interstitial site captures 
three vortices in a trimer state and all, rather than half, 
of the small interstitial sites capture one vortex. This 
state was termed the "second" matching field in Ref. [3l[ . 
Since all of the small interstitial sites are now occupied, 
the sixfold symmetry of the potential in the large intersti- 
tial sites is restored, resulting in two degenerate orienta- 
tions for each trimer. The interaction between neighbor- 



ing trimers is not strong enough to induce a global orien- 
tation of the trimers, so as shown in Fig. nUTa.b) there is 
no long range trimer orientational order at B/B^ = 8/3, 
just as there was no long range orientational order at 
B/Bf = 2.0 in Fig. ifod). There is also only a small 
peak in f c /f p at B/Bf = 8/3 as seen in Fig. UJ^b). 

For B/Bf = 3.0, illustrated in Fig. HDTcd). the large 
interstitial region captures four vortices which form a 
quadrimer state. The remaining interstitial vortices fill 
all of the small interstitial sites. The quadrimers are ori- 
cntationally ordered and have an effective herringbone 
structure in which the quadrimers are tilted in opposite 
directions from one row to the next. This state might 
more appropriately be called a double herringbone, since 
the herringbone state is composed of single dimers, while 
here two dimers have been paired to form a quadrimer. 
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As mentioned earlier, herringbone states appear in sur- 
face physics for the deposition of molecular dimers on 
triangular substrates [35j | . We are not aware of any ob- 
servation of the double herringbone structure for surface 
ordering. Such states might occur for the deposition of 
cubic molecules on triangular surfaces. Fig. [7] shows that 
there is a prominent peak in f c /f p at B/B^ — 3.0 due to 
the overall orientational ordering of the quadrimer state. 

An ordered state also forms at B/B^ = 10/3 as shown 
in Fig. fTUf e.f). Each large interstitial region captures five 
vortices which form a pentagon and the smaller intersti- 
tial sites all capture one vortex. The interstitial pen- 
tagons are aligned in a single direction, as indicated in 
Fig. [He). In surface physics it is uncommon to consider 
fivefold symmetric molecules that form a pentagon; how- 
ever, if such ring type molecules occur and lie flat on an 
atomic surface, then a similar ordering may appear. 

There is no appreciable peak in f c /f P at B/B^ = 11/3 
in Fig. [7J The vortex configuration corresponding to this 
field is plotted in Fig. fTTT a.b). All of the small interstitial 
sites are occupied while some of the large interstitial sites 
capture seven interstitial vortices and the rest capture 
five vortices. The resulting vortex structure has no long 
range orientational order. In general we do not observe 
states where six vortices are captured in the large intersti- 
tial regions. The triangular six- vortex configuration that 
appears in the honeycomb pinning lattice in Fig. Etc) at 
B/B^ — 4 is not stable for the kagome pinning lattice 
due to the presence of the vortices in the small inter- 
stitial sites. At B/B£ = 4.0, which was termed the 
"third" matching field in Ref. [Hj], Fig- [ffic,d) shows 
that each large interstitial region captures seven vortices 
and the remaining interstitial vortices fill the small inter- 
stitial sites. The overall vortex lattice is triangular and a 
prominent peak in f c / f p is observed at this filling as seen 
in Fig. [7J For higher fillings we find additional ordered 
and disordered states and in Fig. fTTT e.f) we show the case 
of B/B^ = 16/3 where a triangular vortex lattice forms. 
At very high fillings B/B^ 3> 1, depinning occurs via a 
shearing motion of the interstitial vortices. 

These results indicate that the kagome pinning ar- 
ray produces peaks in the critical current at most fields 
B/B^ = n/3 for n > 3. In some cases such as 
B/B^ =8/3 and B/B^ = 11/3, the peaks are missing 
or strongly reduced due to a lack of orientational order- 
ing in the vortex configuration. The general behavior for 
both the honeycomb and kagome pinning arrays is very 
similar in that n-mer states with various types of orien- 
tational ordering form in the large interstitial regions of 
the pinning lattice. 



V. VORTEX PLASTIC CRYSTAL STATES 

As noted in Section III(C), the orientationally ordered 
dimer and higher order n-mer states that form in the 
large interstitial regions of the honeycomb and kagome 
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FIG. 11: Left panels: Vortex positions (black dots) and pin- 
ning site locations (open circles) in a portion of the sample 
for the kagome array from Fig. [7] Right panels: Vortex po- 



sitions only. (a,b) B/B$ = 11/3. (c,d) B/B 



(e,f) 



B/B$ = 16/3. 



pinning lattices are very similar to the colloidal molecular 
crystals studied recently [25J, |26|, |28| . Colloidal molecular 
crystals undergo a thermal disordering transition to a 
state in which the orientational ordering of the n-mers is 
lost when the n-mers remain localized at the lattice sites 
but begin to rotate freely [ID, [26|, [28| . Here we show that 
the vortex n-mers exhibit a similar rotational melting 
behavior. 

In Fig. [T27 a) we plot the vortex positions and pin- 
ning site locations for the honeycomb pinning array in 
Fig.[3a,b) at B/Bf = 2 and temperature F T = 1.56. At 
F T = in Fig. a), the dimers are all aligned in a single 
direction, but at F T = 1.56 in Fig. [T27 a) the dimer align- 
ment is lost, although the dimers remain confined to the 
interstitial sites. The vortex trajectories, illustrated in 
Fig. [T2T b) for a period of 10 3 simulation time steps, indi- 
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FIG. 12: (a) The vortex positions (black dots) and pinning 
site locations (open circles) in a portion of the sample for the 
honeycomb pinning array in Fig. GJa,b) with f p = 0.5/o and 
B/Bf = 2.0 at a temperature of F T = 1.56. (b) The same as 
in (a) with vortices (white dots), pinning sites (open circles), 
and vortex trajectories (black lines) over a fixed time period 
of 10 3 simulation time steps. 



cate that while the pinned vortices move a small amount, 
each dimer is undergoing rotational motion between the 
three low energy orientations. We find that there is some 
correlation in the dimer motions; however, over long dis- 
tances the true long range orientational ordering is lost. 
We refer to this finite temperature state with no orienta- 
tional order as a "vortex plastic crystal" since it is sim- 
ilar to the colloidal plastic crystal phases. In the vortex 
dimer case, there are two species of vortices: the intersti- 
tial dimers and the single vortices trapped at the pinning 
sites. In the colloidal system of Refs. [25|, [2(| [28[, the 
egg-carton substrate does not allow interstitial colloids 
to exist so there is only a single species of n-mer states. 

In Fig. [T37 a) we plot the fraction of orientationally or- 
dered dimers P as a function of temperature F T for the 
honeycomb pinning lattice at B/Bf = 2 from Fig. QH 
To determine P, for each dimer we identify the angle 
a which the line connecting the two vortices compos- 
ing the dimer makes with the positive x axis, where 
< a < 180°. We then assign a state di to each dimer 
of 1, 2, or 3 depending on whether a is closest to 30°, 
90°, or 120°. The ordered fraction P is given by P = 

■di)) 
For 



max(£ti 6(1 *)> *(? ~ * 



,£ti<5(3- 

where is the number of dimers in the sample. 
F T < 1.25 the dimers remain orientationally ordered 
with P = 1 and form the vortex molecular crystal state 
shown in Fig. Efa). For 1.25 < F T < 2.3, the orienta- 
tional order is thermally destroyed and the system enters 
the vortex plastic crystal state. For F T > 2.3 the dimers 
break apart, P is no longer defined, and there is vortex 
diffusion throughout the entire system. 

Figure [K^b) shows the depinning force f c /f P as a 
function of F T for the honeycomb pinning lattice at 
B/Bfp = 2 from Fig. UMa,). The orientationally or- 
dered dimer state has a well defined depinning threshold 
which decreases monotonically with increasing F T . Once 
the dimer orientational order is lost, f c /f P undergoes a 




FIG. 13: (a) Fraction P of orientationally ordered dimers 
versus temperature F T for the honeycomb pinning lattice in 
Fig.EK b ) with f p = 0.5/o and B/Bf = 2.0. (b) The depin- 
ning force fc/fp vs F T for the same sample. 



corresponding pronounced drop. This is consistent with 
the results presented earlier in which vortex states that 
lack orientational order have a significantly lower value 
of fc/fp than orientationally ordered states. We note 
that in the experimental studies on honeycomb arrays of 
Ref. [33j], as the temperature was increased the match- 
ing effects at B/B^ = 3.5 and 4.5 were lost while the 
matching effects at B/B^ = 1 and 2 persisted. This is 
consistent with our finding that a significant drop in the 
strength of the matching effect appears at the melting 
of the vortex molecular crystals, which we expect to be 



present at B/B¥ = 3.5 and 4.5. 



VI. 



VORTEX MOLECULAR CRYSTAL PHASE 
DIAGRAM 



A. Honeycomb Dimer State at B/B. 



By conducting a series of simulations and measuring P 
and the diffusion D we determine the phase diagram of 
the different phases for the B/B^ = 2.0 dimer state of 
the honeycomb pinning array. The diffusion is given by 



D = 



\R t (t + dt) -Rj(t)| 
dt 



(4) 
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FIG. 14: The temperature F T vs pinning strength f p phase 
diagram for the honeycomb pinning array at B/B^ = 2.0 in 
the dimer state. The four phases include: the vortex molec- 
ular crystal phase illustrated in Fig. [SJa,b), a vortex plastic 
crystal phase shown in Fig. I13f a,b), a modulated vortex liq- 
uid phase where motion occurs throughout the sample, and a 
partially pinned triangular vortex lattice that forms at low f p 
and low temperature which is described in Fig. [15] Circles: 
Onset of significant diffusion as determined from the diffu- 
sion measurement D. Squares: Loss of orientational order 
as determined from P, the fraction of orientationally ordered 
dimers. 



with dt = 1000 simulation time steps. The resulting 
phase diagram as a function of temperature F T and pin- 
ning strength f p is given in Fig. [14] The behavior of the 
sample at f p — 0.5/o was shown in Figs . [P2l and 1 1 31 where 
a vortex molecular crystal forms for < F T < 1.25, a 
vortex plastic crystal is present for 1.25 < F T < 2.3, and 
at F T > 2.3 the dimer states break apart and the vor- 
tices diffuse throughout the sample in a modulated liquid 
state induced by the substrate. Figure [T4l indicates that 
there is no vortex plastic crystal phase for f p < 0.35/o 
but that the vortex molecular crystal melts directly to 
the modulated vortex liquid for 0.15/o < f p < 0.35/o- 

For f p < 0.15/o, the pinning force is not strong enough 
to overcome the elastic energy of the vortex lattice and 
the low temperature ground state is not the dimer state 
shown in Fig. [5fa,b) but a partially pinned triangular 
vortex lattice. This state is illustrated in Fig. [TST a.b) 
a t fp = 0.15/o and F T = 0. The triangular vortex lat- 
tice still shows a partial commensuration effect with the 
substrate and is aligned in such a way that half of the 
pinning sites are occupied by a vortex. The unpinned 
vortices form a kagome structure oriented 30° from the 
original triangular lattice, where here the pinned vortices 
take the place of the missing pinning sites in the kagome 
lattice. 

The melting temperature in Fig. [H] from the partially 
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(a) x (b) x 

FIG. 15: (a) Vortex positions (black dots) and pinning site 
locations (open circles) in a portion of the sample for the 
honeycomb pinning array in Fig.[l4]at B/B^ — 2.0 and f p = 
0.15/o in the partially pinned triangular lattice phase at F T = 
0. (b) Vortex positions only. 

pinned triangular lattice to the modulated vortex liquid 
increases with increasing f p . The existence of pinned 
vortices in the triangular vortex lattice effectively stiff- 
ens the lattice, and this effect becomes more pronounced 
as f p increases. The transition temperature from the vor- 
tex plastic crystal state to the modulated vortex liquid 
also increases with increasing f p since a higher tempera- 
ture is required to enable the vortices to escape from the 
stronger pinning sites. 

An interesting feature in the phase diagram of Fig. [14] 
is the fact that the transition temperature between the 
vortex molecular crystal and vortex plastic crystal states 
decreases with increasing f p . The same trend was ob- 
served for colloidal molecular crystals as a function of 
temperature versus substrate strength [25|, [U UK]. In 
the colloidal case, each dimer is trapped in a substrate 
minima, and as the substrate strength increases, the dis- 
tance between the two colloids in each dimer decreases. 
As a result, the strength of the effective quadrupole inter- 
action between the dimers decreases, lowering the tran- 
sition temperature. For the case considered here, the 
vortex dimers in the honeycomb pinning lattice do not 
sit in the pinning sites and are not directly affected by 
the increase of f p . The dimers do, however, interact with 
the vortices which are trapped by the pinning sites, and 
these pinned vortices become less mobile and less able to 
fluctuate as f p increases. At low f p , the pinned vortices 
that are closest to the ends of each dimer move to the 
outer edges of the pinning sites in order to sit as far as 
possible from that dimer. This produces a natter inter- 
stitial confining potential along the length of the dimer in 
the center of the honeycomb plaquette and permits the 
interstitial vortices that form the dimers to move further 
away from each other. For higher f p , the pinned vor- 
tices are shifted toward the center of the pinning sites 
and closer to the dimer. This produces a stronger confin- 
ing force on the interstitial dimer and brings the vortices 
that form the dimer closer together, reducing the effective 
quadrupole interaction between neighboring dimers and 
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FIG. 16: (a) The vortex positions (black dots) and pinning 
site locations (open circles) in a portion of the sample for the 
honeycomb pinning array in Fig. Oc,d) with f p = 0.5/o and 
B/Bf — 2.5 at a temperature of F T = 1.1. (b) The same as 
in (a) with vortices (white dots), pinning sites (open circles), 
and vortex trajectories (black lines) over a fixed time period 
of 10 3 simulation time steps. 



lowering the transition temperature between the vortex 
molecular crystal and vortex plastic crystal states. This 
transition line is likely to saturate at very high f p when 
the pinned vortices are constrained to sit at the very cen- 
ter of the pinning sites and cannot adjust their positions 
in response to the orientation of the neighboring dimers. 
We note that it is possible that for very high values of 
f p , a new interstitial liquid phase could form in which 
the dimers break apart and the interstitial vortices hop 
from one interstitial region to another while the vortices 
in the pinning sites remain immobile. 

It is beyond the scope of this work to determine the 
exact nature of the transitions in the phase diagram of 
Fig. [Mj however, since the overall system appears to be 
very similar the colloidal case, we can argue from the 
results for the colloidal system that the vortex molecu- 
lar crystal to vortex plastic crystal transition is probably 
Ising-like. The vortex plastic crystal to modulated vor- 
tex liquid state transition is mostly likely an activated 
crossover. 



B. Honeycomb Trimer State at B / B" = 2.5 

In order to determine how general the features of the 
phase diagram in Fig. [14] are, we consider the case of 
the honeycomb lattice at a field of B/B^ — 2.5 where 
an ordered arrangement of trimers occurs as seen in 
Fig.[5Jc,d). We find the same general melting behavior as 
for the dimers. In Fig.fTWa) we illustrate the vortex plas- 
tic crystal at F T = 1.1. Here each pinning site captures 
one vortex and the interstitial trimers do not have ori- 
entational ordering. In Fig. RTOf b) the vortex trajectories 
indicate that the trimers are rotating in a manner similar 
to the dimers in Fig. [T2] For higher temperature F T , the 
trimers break apart and we observe diffusion throughout 
the entire sample. 




Polycrystal 



Vortex Molecular Crystal 
i I i 
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FIG. 17: The temperature F T vs pinning strength f p phase 
diagram for the honeycomb pinning array at B/B^ = 2.5 
in the trimer state. The vortex molecular crystal phase is 
illustrated in Fig. [SJc,d). The vortex plastic crystal phase is 
shown in Fig. 1161 In the modulated vortex liquid, there is 
diffusion throughout the entire sample, while at low f p and 
F T a partially pinned poly crystalline triangular lattice forms 
which is described in Fig. 1181 Circles: Onset of significant 
diffusion. Squares: Loss of trimer orientational order. 



By conducting a series of simulations for varied tem- 
perature F T and pinning strength f p and measuring the 
fraction of aligned trimers and the diffusion, we map out 
the phase diagram for the trimer state at B/B^ = 2.5 
as shown in Fig. |T7j The general features of the phase 
diagram at B/B^ = 2.5 are similar to the phase dia- 
gram at B/B^ = 2.0 in Fig. [HI however, there are some 
noticeable differences. For a given pinning strength, the 
trimers disorder at a significantly lower temperature than 
the dimers. For example, at f p = 0.5/o, the trimers lose 
orientational order at F T = 0.5 while for dimers the ori- 
entational order persists up to F T — 1.25. This reflects 
the fact that the effective multipole interaction between 
trimers is weaker than the effective quadrupole interac- 
tion between dimers. We also note that at F T = 0.0 
we observe a peak in f c /f p at B/B^ = 2.5, shown in 
Fig. [21 while the experiments of Ref. [33[ did not find a 
peak at this filling. This could be due to the fact that 
the trimer melting temperature is relatively low, and the 
experiments may have been performed above the trimer 
melting temperature. 

For f p < 0.2/o and low temperatures, the ordered 
trimer state is lost and the vortices form a partially or- 
dered triangular lattice as shown in Fig. [TH] At B/B^ = 
2.5, the vortices are unable to simultaneously sit in a 
triangular lattice and align with the pinning sites, as 
in the partially pinned triangular lattice of Fig. [15] at 
B/Bjp = 2.0. Instead we find a polycrystalline state 
composed of a triangular lattice that contains disloca- 
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(a) 



FIG. 18: (a) Vortex positions (black dots) and pinning site 
locations (open circles) in a portion of the sample for the 
honeycomb pinning array in Fig. [17] at B/B^ — 2.5 and f p — 
O.l/o in the polycrystalline vortex state at F T = 0. (b) Vortex 
positions only. 



tions and grain boundaries. The melting line between 
the polycrystalline state and the modulated vortex liq- 
uid in Fig. [17] at B/B^ = 2.5 decreases in temperature 
with increasing f p . This is the opposite of the behavior of 
the melting line between the partially pinned triangular 
lattice and the modulated vortex liquid at B/B^ = 2.0 
shown in Fig. 1 141 and occurs because the increasing f p 
leads to an increase of the polydispersity in the lattice, 
depressing the melting temperature. At f p — 0.0 and 
B/B^ = 2.5, a dislocation-free triangular lattice forms 
which has a higher melting temperature than the defected 
lattice that appears at finite pinning strength. 

We have also studied the effect of temperature on the 
orientational ordering of the n-mer states for B/B^ > 3 
in the honeycomb lattice at fixed f p = 0.5 (not shown). 
In general we find that at fields with no long-range ori- 
entational order, such as B/B^ ~ 3, the n-mers undergo 
thermally induced rotations at any finite temperature so 
there is only a vortex plastic crystal phase which melts 
to the modulated vortex liquid. At the other orienta- 
tionally ordered fillings, we observe a finite temperature 
transition from an ordered vortex molecular crystal state 
to a vortex plastic crystal state. 

We have performed finite temperature simulations for 
the kagome system as well and find the same general 
results (not shown) as in the honeycomb system. At ori- 
entationally ordered fillings, there is a low temperature 
ordered vortex molecular crystal state which melts into 
an intermediate vortex plastic crystal state. Additionally, 
the orientationally disordered states at B/B^ = 8/3 and 
11/3 show only vortex plastic crystal and vortex modu- 
lated liquid phases at finite temperature. 

Our results agree with the previous study of thermal 
motion of vortices on kagome pinning arrays. In Ref. [3l| . 
at the second matching field B/B^ = 2, which corre- 
sponds to our B/B^ = 8/3, the interstitial vortex tri- 
angles did not form an orientationally ordered state, but 
instead created what was termed a kagome state at fi- 



nite temperatures. This kagome state corresponds to 
our general class of vortex plastic crystals in which the 
n-mers are orientationally disordered. We note that at 
B/B^ = 8/3, we observe no large peak in f c /f p as shown 
in Fig. [7] and the interstitial vortex trimers have no long- 
range orientational order as seen in Fig. UOf a). In con- 
trast, at B/B£ = 7/3 we find an orientationally ordered 
vortex molecular crystal state in which all of the trimers 
are aligned, shown in Fig. [He), which has a finite tem- 
perature transition to a vortex plastic crystal. 

The vortex plastic crystals have similarities to the re- 
cently proposed vortex Peierls states for vortices in a 
Bose-Einstein condensate interacting with a co-rotating 
periodic optical lattice [23[ . The optical lattice structure 
is in fact a kagome array; however, in Ref. [23| there is an 
additional potential minima imposed on the center of the 
kagome plaquettes. The vortex trimer state which forms 
at a filling of 1/3 on the dual dice lattice tunnels between 
the two degenerate configurations. This state is very sim- 
ilar to the vortex trimer state we observe in which the 
trimers are undergoing thermally induced rotations, as 
shown in Fig. 1161 



VII. EFFECT OF PINNING STRENGTH 

The strength f p of the pinning sites determines 
whether it is possible for a vortex molecular crystal state 
to form. In Sec. VI we showed that in the honeycomb pin- 
ning lattice at B/B^ = 2.0, low f p and low temperature, 
the partially pinned triangular vortex lattice illustrated 
in Fig. [15] appears instead of the vortex dimer molecular 
crystal shown in Fig. (Ha). This results when the elas- 
tic energy of the vortex lattice, which favors a triangular 
vortex configuration, overcomes the energy of the pin- 
ning sites. Similarly, at B/B^ = 2.5 the partially pinned 
polycrystalline vortex lattice seen in Fig. [reforms instead 
of a vortex trimer molecular crystal at low f p . In general, 
partially pinned phases will not occur in superconductors 
with arrays of holes because the pinning is too strong, 
but may form for the weaker pinning found for vortices 
in superconductors with blind hole arrays or for colloids 
on optical trap arrays. Previous work has shown that a 
competition between the symmetry of the pinning lattice 
and a triangular lattice of interacting particles can lead 
to structural transitions of the particle lattice. Transi- 
tions from square to partially pinned vortex lattices have 
been observed in vortex simulations [2l|, [37| , while par- 
tially pinned phases have recently been demonstrated ex- 
perimentally for macroscopic Wigner crystals in square 
pinning arrays [30] . As the pinning strength is reduced at 
the first matching filling, Ref. [30( shows that a transition 
occurs from a square Wigner lattice where each pinning 
site captures one charge to a partially pinned triangular 
lattice. A similar transition appears at the second match- 
ing filling. If the substrate had been triangular, then at 
the first matching filling there would have been no struc- 
tural transition of the charge lattice since its symmetry 
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FIG. 19: f c vs B/Bf for the honeycomb pinning array from 
Fig.CHwith f p = 0.35/o, 0.175/o, and 0.1/ , from top to bot- 
tom. The results for f p — 0.5/o appear in Fig. [2] 



would coincide with that of the pinning lattice. At the 
second matching filling on a triangular array, the parti- 
cles form a honeycomb lattice when the pinning is strong 
0, but as the pinning strength is weakened, a transition 
to a partially pinned phase should occur that allows the 
particle lattice to have a more triangular ordering. For 
vortices in honeycomb and kagome pinning arrays, we 
expect a transition from a completely pinned phase to a 
partially pinned phase to occur as a function of pinning 
strength even at the first matching field. 

Since there are a rather large number of possible dif- 
ferent vortex configurations shown in Sec. Ill and Sec. 
IV, we focus here on three cases which indicate the gen- 
eral behavior for both the honeycomb and kagome arrays 
at the different matching fields as the pinning strength 
is varied. In Fig. [TH] we plot f c vs B/B^ at pinning 
strengths of f p = 0.1/o, 0.175/o, and 0.35/o f° r the 
honeycomb pinning array from Fig. [21 which contains 
the results at f p = 0.5/o- The vortex configurations at 
f p = 0.35/o are the same as those at f p — 0.5/o described 
in Section III, and the depinning force f c at B/B^ =1/2 
and 1.0 is higher than f c at B / B" = 1.5 and 2.0. As f p 
is reduced below 0.35/o, the overall f c at all the fields 
decreases, but not uniformly, as f c at B/B^jf = 1.0 drops 
below the value of f c at B/B^ = 0.5 and 1.5. This 
crossover in f c with decreasing pinning force results when 
the vortex lattice at B/B^ = 1 undergoes a transition 
from a fully pinned honeycomb lattice such as that il- 
lustrated in Fig. [3|c,d) to the partially pinned distorted 
square lattice shown in Fig. I2"07a.b) at f p = 0.175/o- 
Here, although there are an equal number of pins and 
vortices, half of the pinning sites are occupied while the 
other half are empty, resulting in a significant decrease 
of the depinning threshold. The elastic energy of the dis- 
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FIG. 20: (a) The vortex positions (black dots) and pinning 
site locations (open circles) in a portion of the sample for 
the honeycomb pinning array in Fig. \T§\ at B/B^ — 1.0 and 
f p = 0.175/o where a distorted square vortex lattice forms, 
(b) Vortex positions only. For higher f p , all the vortices are 
located at pinning sites and the state shown in Fig. HKc,d) 
occurs. 
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FIG. 21: (a) The depinning force f c vs f p for the honeycomb 
pinning array at B/B^f — 1/2 (squares), B/B^f = 1.0 (cir- 
cles), and B/B^ — 1.5 (triangles). At f p < 0.35, the vortex 
lattice at B/Bf = 1 is a distorted square lattice rather than a 
honeycomb lattice, and therefore f c for B/B^ — 1 drops be- 
low f c for B/B^ = 1/2. (b) The same as in (a) for B/Bf = 1 
(circles) and B/B^ — 1.5 (triangles) but extended to higher 
f p , showing that there is a saturation in / c for B/B^ = 1.5. 



torted square vortex lattice at f p < 0.35/o is lower than 
the honeycomb vortex lattice that forms at higher f p . 

In Fig. l2"TT a) we plot f c vs f p for the honeycomb ar- 
ray at B/Bf = 1/2, 1.0, and 1.5. At B/B^ = 1/2, the 
triangular vortex lattice illustrated in Fig. 0Ja,b) forms. 
There is no transition to a partially pinned state as f p de- 
creases at this field since there is no competition between 
the vortex lattice symmetry and the pinning symmetry. 
All the vortices remain trapped at pinning sites for ar- 
bitrarily low f p . Since the vortex lattice is completely 
symmetrical the depinning force f c is directly propor- 
tional to f p , as shown in Fig. [231a). At B/B^ = 1.0, 
all the vortices are pinned at the pinning sites in a hon- 
eycomb lattice for f p > 0.35/o and f c increases linearly 
with f p . For f p < 0.35/o the system enters the partially 
pinned phase illustrated in Fig. [201 At the same time, 
f c drops abruptly to a lower value and then decreases 



17 



more slowly with decreasing f p . The vortex lattice is al- 
ways triangular at B/B^ = 1.5, as seen in Fig. Efe,f), 
and there is no sharp change in f c vs f p since there is no 
structural transition in the vortex lattice when f p is var- 
ied. For f p < 0.35, f c is higher at B/B^ = 1.5 than at 
B/B^ = 1.0, as also indicated in Fig. [113 The depinning 
force at B/B% = 1.5 increases more slowly than linearly 
with f p , which is more clearly seen in Fig. [2lT b) where 
the range of f p is extended to 2.0/ . Here, the depin- 
ning force is determined by the 1/3 of the vortices which 
are not confined at pinning sites. As f p increases, the 
vortices at the pinning sites are more strongly pinned; 
however, the caging potential which pins the interstitial 
vortices saturates. 

At B/B^ = 2.0 we have already shown in Section VI 
(B) that a transition to the partially pinned state illus- 
trated in Fig. [15] occurs as f p decreases. As a result, 
there is a jump in f c near /„ = 0.35/o (not shown) simi- 
lar to what we find for B/B^ = 1 in Fig. 121( a). For high 

f pi the f c versus f p curve for B/B^ = 2 resembles that 
at B/B^ = 1.5 shown in Fig. l2~lTb), with a saturation 
in the depinning force for the interstitial vortex dimers. 
We thus expect two dominant behaviors of f c versus f p 
at higher matching fields. For those fields where the 
overall vortex lattice is triangular, such as B/B^ = 1.5 
and 4.5 in the honeycomb pinning array (Fig. Hfe,f) and 
Fig. E](e,f), respectively) as well as B/B^ = 4/3, 4, and 
16/3 in the kagome array (Fig. [5](e,f), Fig. [TlTc.d). and 
Fig. [TlTe.f). respectively), f c saturates at high values of 
f p but has no transition at low values of f p . For the 
matching fields with nontriangular vortex configurations, 
a transition to a more triangular state at low f p occurs 
which is accompanied by a drop in f c , and in addition 
there is a saturation of f c at high f p which always oc- 
curs in the presence of interstitial vortices. It is also 
possible that there could be more than one structural 
transition as f p decreases. For example, we have shown a 
transition from a honeycomb vortex lattice to a distorted 
square lattice (Fig. |2"0)) in the honeycomb pinning lattice 
at B/B^ = 1. At extremely low f p , a second transition 
to a completely triangular vortex lattice should occur. 



VIII. DISCUSSION 

In this work we only considered the case where each 
pinning site could capture at most one vortex; however, 
other types of vortex phases and commensurability ef- 
fects may occur in real samples if multiple vortices are 
captured by the pins. If the pinning is very strong, the 
first few matching fields would correspond to the pres- 
ence of multiquanta vortices at the pinning sites and no 
interstitial vortices. The vortices would start to enter the 
interstitial regions only once the pinning site saturation 
field is reached. For kagome and honeycomb lattices of 
this type, strong commensurability effects would appear 
at only the integer matching fields below the pinning sat- 



uration field, and would cross over to matching effects at 
the n/2 or n/3 fields once the pinning saturation occurs. 
Alternatively, individually quantized vortices may form 
at low fields, while multiply quantized vortices may ap- 
pear in the pinning sites only at high fields where the 
vortex lattice constant is small. In such a case, noninte- 
ger matching effects would appear at low fields followed 
by a crossover to integer matching effects at higher fields. 
In Ref. 0, where continuum simulations of square pin- 
ning arrays were performed, it was proposed that multi- 
ple interstitial vortices can merge to form a single giant 
interstitial vortex. In the case of the honeycomb and 
kagome arrays, singly quantized vortices could occupy 
the pinning sites while multiply quantized vortices sit in 
the large interstitial sites. In this scenario the overall 
vortex lattice would be triangular and commensurability 
peaks would be observable at all fields n/3 or n/2. If mul- 
tiple quantization of the interstitial vortices occurs, there 
would be no finite temperature phase transition from a 
vortex molecular crystal to a vortex plastic crystal and 
all of the matching peaks would vanish at the same rate 
with temperature. Additionally, since the vortex lattice 
symmetry would be triangular at every matching field, 
there would be no missing matching peaks due to the 
formation of vortex plastic crystals. These effects are 
not seen in the experiment of Ref. [33| . suggesting that 
multiple quantization of vortices is not occurring. 

As stated previously, many of the orientationally or- 
dered n-mer states are related to various types of Ising 
and three-state Potts models. It may be possible to cre- 
ate honeycomb and kagome pinning arrays that have an 
additional anisotropy in one direction. This would bias 
some of the degenerate directions of the n-mer states so 
that only certain directions of the n-mer ordering would 
occur. Also, highly frustrated states could be created 
such as incommensurate fillings composed of mixtures of 
two different n-mer species. Ordered colloidal molecular 
crystal alloy states formed by mixtures of different n-mers 
have been proposed in Ref. (38j . and similar states may 
occur for the vortex system. It is likely that such mix- 
tures would have extremely long relaxation times to reach 
the ordered ground state. These mixtures might produce 
glassy dynamics and have interesting time-dependent or 
history-dependent properties. 

In addition to allowing the creation of vortex molecular 
crystal and vortex plastic crystal states, honeycomb and 
kagome pinning arrays may also be valuable for the gen- 
eral enhancement of pinning. If pinning sites are placed 
in a honeycomb or kagome arrangement, higher pinning 
can be achieved at matching fields above the first match- 
ing field compared to a triangular pinning array with 
an equal number of pinning sites. This is illustrated in 
Fig. |2"21 where f c / f p is plotted as a function of matching 
field for a triangular, kagome, and honeycomb pinning 
array of equal pinning strength f p = 0.5/o- Below the 
first matching field there is little difference between the 
depinning force for the three arrays; however, above the 
first matching field the honeycomb array shows a pro- 
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FIG. 22: (Color online) f c /fp versus B/B^ for the triangular 
pinning array (red medium line) , versus B/B$ for theka gome 
array (dark line), and versus B/B^ for the honeycomb array 
(light line) at f p = 0.5/o. 



nounced enhancement of f c over the triangular array and 
the kagome array also shows a smaller enhancement of f c 
compared to the triangular array. This suggests that for 
different arrays with the same number of pinning sites, 
the honeycomb and kagome pinning arrays have an over- 
all higher critical current than triangular pinning arrays. 



IX. CONCLUSION 

We have used numerical simulations to study the vor- 
tex states in honeycomb and kagome pinning arrays in 
the limit where only one vortex is captured at each pin- 
ning site. For the honeycomb arrays, we find pronounced 
matching effects in the form of peaks in the depinning 
force at most fields B/B^ = n/2 where n > 2, while 
for kagome arrays pronounced matching effects occur at 
most fields B/B^ = n/3 where n > 3. This is in con- 
trast to the purely triangular pinning arrays which have 
prominent matching effects only at B /B^ = n or at fields 
below the first matching field. For the honeycomb pin- 
ning array, a variety of novel vortex molecular crystal 
states occur. Here, multiple interstitial vortices occupy 
the large interstitial regions of the honeycomb pinning 



lattice and form effective dimer, trimer and higher order 
n-mer vortex states. The n-mers interact via an effective 
quadrupole or higher order pole moment which can lead 
to an overall n-mer orientational ordering. This type of 
?i-mcr ordering is similar to the recently studied colloidal 
molecular crystal states observed for repulsively inter- 
acting colloidal particles on periodic substrates. For the 
honeycomb array, at some of the matching fields the ri- 
mers do not order and there is no prominent peak in the 
depinning force at these fields. Our results agree well 
with recent experiments on honeycomb pinning arrays 
where strong peaks are observed at B/B^ = n/2 fillings 
with missing or weak peaks for the fields at which we ob- 
serve disordered n-mer configurations. For the kagome 
pinning arrays we observe similar dimer, trimer and 
n-mer ordered states which form herringbone or other 
ferromagnetic-like configurations. Most of the missing 
commensuration peaks in the B/B^ — n/3 sequence 
correspond to the lack of n-mer ordering at zero temper- 
ature. We term the orientationally ordered vortex n-mer 
states "vortex molecular crystals." At finite temperature 
we demonstrate the existence of a transition from an ori- 
entationally ordered vortex molecular crystal state to a 
state where the n-mers are rotating and lose their relative 
orientational ordering. We refer to this disordered state 
as a vortex plastic crystal. At higher temperatures there 
is a crossover to a modulated vortex liquid state where 
there is diffusion throughout the entire system. When 
the pinning strength is weak, we find that the vortex 
molecular crystal phases undergo a transition to partially 
pinned phases where the vortex lattice has triangular or 
partially triangular ordering and only a portion of the 
pinning sites are occupied. These results suggest that the 
vortex molecular crystal states have many similarities to 
the recently studied colloidal molecular crystal states and 
that the vortex molecular crystal to vortex plastic crys- 
tal transitions can be mapped to various types of spin 
systems such as Ising and Potts models. 
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